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ABSTRACT 



We present new HI observations of the nearby dwarf galaxy NGC 3741. This galaxy 
^ I has an extremely extended HI disk, which allows us to trace the rotation curve out to 

unprecedented distances in terms of the optical disk: we reach 42 B-band exponential 
scale lengths or about 7 kpc. The HI disk is strongly warped, but the warp is very 
symmetric. The distribution and kinematics are accurately derived by building model 
data cubes, which closely reproduce the observations. In order to account for the 
\^ I observed features in the data cube, radial motions of the order of 5. . .13 km s^^ are 

. needed. They are consistent with an inner bar of several hundreds of pc and accretion 

of material in the outer regions. 

The observed rotation curve was decomposed into its stellar, gaseous and dark 
components. The Burkert dark halo (with a central constant density core) provides 
O ' very good fits. The dark halo density distribution predicted by the ACDM theory fails 

to fit the data, unless NGC 3741 is a 2.5-a exception to the predicted relation between 
^2 ' concentration parameter and virial mass and has at the same time a high value of the 

^ ' virial mass (though poorly constrained) of 10^^ M©. Noticeably, MOND seems to be 

consistent with the observed rotation curve. Scaling up the contribution of the gaseous 
■ disk also gives a good fit. 

5h ' Key words: galaxies: kinematics and dynamics - galaxies: spiral - dark matter 



1 INTRODUCTION 

The measurement of rotation curves of disk galaxies has 
been one of the most powerful tools to study dark matter, 
its content relative to baryons, and its distribution. In par- 
ticular, dwarf galaxies are good candidates for dark-matter 
studies as their kinematics is generally dominated by dark 
matter down to small galactocentric radii (Persic, Salucci & 
Stel 1996), thus providing an almost clean measurement of 
the dark-matter contribution to the observed rotation curve 
and hence of its density profile. Therefore, they represent 
an important test for cosmological models, because numeri- 
cal simulations performed in the framework of the currently 
favoured theory of structure formation in the Universe, A 
Cold Dark Matter (ACDM), predict a well-defined density 
profile for dark matter in virialised structures, the NFW 
profile (Navarro, Frenk and White 1996). Subsequent mod- 
ifications and refinements via improved numerical resolu- 
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tion (Moore et al. 1999, Navarro et al. 2004) do not change 
one of the main predicted features: a central density cusp, 
p(r) oc for r ^ 0, with 1.0 < a < 1.5. 

In the recent past, several studies were devoted to the 
analysis of rotation curves of dwarf galaxies (e.g. van den 
Bosch & Swaters 2001, Borriello & Salucci 2001, Swaters 
et al. 2003, Weldrake, de Blok & Waher 2003, Simon et 
al. 2005, Gentile et al. 2005), in particular to test whether 
the density profile predicted by ACDM models is actually 
observed. Even though the debate as to whether observa- 
tions can put strong constraints on the dark matter density 
profiles is still ongoing, the results are that in most cases 
the observations prefer halos with a constant-density core, 
p(r) ~ po for r 0. A possible explanation to reconcile the 
observations and predictions involves non-circular motions 
induced by halo triaxiality (Hayashi et al. 2004, Hayashi & 
Navarro 2006); the authors show the effects of non-circular 
motions on the observed minor axis kinematics for some 
viewing angles. Gentile et al. (2005) analysed the velocity 
field of DDO 47 (of which the rotation curve is best fitted 
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Figure 1. HI total intensity map from the liigli-resolution data 
cube. Contours are (3, 4.5, 6.75, ...) X lO^" atom cm . The beam 
is shown in the bottom left corner. 
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Figure 2. HI total intensity map from the low-resolution data 
cube (contours) superimposed onto an optical image (DSS, 
greyscale). Contours are (1.5, 2.25, 3.375, ...) X 10^^ atom cm~-^. 
The radio beam is shown in the bottom left corner. 



by a cored halo) through its harmonic decomposition and 
found non-circular motions of the order of only 2 — 3 km 
s~^, an order of magnitude smaller than the discrepancy 
between the observations and the ACDM predictions. How- 
ever, the effects of triaxiality vary strongly with the viewing 
geometry. A detailed study of the effects of triaxiality with 
all possible viewing angles on the results of harmonic de- 
compositions of velocity fields is yet to come and would be a 
valuable result. It seems however unlikely that the very small 
2-D non-circular motions measured in Gentile et al.(2005) 
may hide an underlying NFW cusp. 

NGC 3741 is a sort of privileged laboratory to investi- 
gate the dark matter distribution in galaxies. It is a nearby 
dwarf irregular galaxy with an absolute blue magnitude of 
-13™13. Its distance (Karachentsev et al. 2004), estimated 
through the tip of the Red-Giant Branch, is 3.0±0.3 Mpc, 
while Georgiev et al. (1997) find a slightly higher value (3.5 
± 0.7 Mpc) from the brightest blue stars; unless stated oth- 
erwise, in this paper we adopt a distance of 3 Mpc, obtaining 
an absolute B-band blue luminosity of 2.7 x 10^ Lq. The 
HI disk of this galaxy was studied by Begum et al. (2005) 
using observations performed with the Giant Metrewave Ra- 
dio Telescope (GMRT), which disclosed the extremely ex- 
tended gaseous disk of NGC 3741. They traced the rotation 
curve out to 38 B-band times the exponential scale length, 
whose value is 10'.' 75 (Bremnes et al. 2000), corresponding 
to 0.16 kpc at a distance of 3 Mpc. We observed this galaxy 
in the HI line with the Westerbork Synthesis Radio Tele- 
scope (WSRT) and reach 42 exponential scale lengths. We 
re-determined the rotation curve by means of the more ef- 
ficient Modified Envelope Tracing method (Gentile et al. 
2004) and data cube modelling. This is to our knowledge 
the most extended rotation curve ever measured in terms of 
the optical size. 



2 OBSERVATIONS AND REDUCTION 

The observations were performed with the WSRT for 12 
hours with the 'maxi-short' configuration. The correlator 
setup was chosen to have 10 MHz total bandwidth di- 
vided into 2048 channels, which were subsequently reduced 
through Hanning smoothing. The final data cubes have a ve- 
locity resolution of 4.1 km s~^. Flagging and calibration were 
performed with the MIRIAD software. A continuum map 
was created from the line-free channels, which was used for 
phase self-calibration. The continuum was subtracted from 
the u-v data using line-free channels and then an HI data 
cube was created, using uniform weighting. The data were 
CLEANed in the regions of the channel maps with HI emis- 
sion and subsequently restored using a Gaussian beam. The 
beam size is 15'.'2 x ll'.'7, which at the distance of NGC 3741 
corresponds to approximately 220 pc x 170 pc. We also cre- 
ated a low-resolution data cube, with a beam of 60" x 60". 
In order to separate real and spurious emission, the data 
cubes were smoothed at about twice the resolution; the re- 
gions in these smoothed cubes with emission lower than 2a 
were blanked. Based on these masked data cubes, total HI 
maps (moment maps) from the high- and low-resolution 
cubes were created: they are shown in Figs. [T] and (2] 

The integrated HI flux is 59.6 Jy km s~^, which corre- 
sponds to an HI mass of 1.3 x 10* Mq, and a total gas mass 
(including primordial He) of Mgas = 1.7x 10* Mq. The value 
of our measured total flux is higher than the single-dish one 
(Schneider et al. 1992), but it is lower than the value quoted 
in Begum et al. (2005). However, it is unlikely that we have 
lost a signiflcant amount of flux due to missing short spac- 
ings, because the largest angular scale imaged by the WSRT 
in the maxi-short configuration is much larger than the size 
of the emission in the channel maps. 
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Figure 3. Parameters of our first attempt to derive tlie rotation 
curve of NGC 3741: the tilted-ring fit to tlie iiigii-resolution ve- 
locity field, with the incUnation fixed and the position angle as a 
free parameter. 



NGC 3741 



o 
o 




11°38™30= 36"15 



35^45= 



35'"30= 



Right Ascension (J2000) 



Figure 4. Comparison between the observed velocity field de- 
rived from the Modified Envelope Tracing method (black con- 
tours, see text) and the modelled one (grey contours), built with 
the output position angles from the tilted-ring fit to the velocity 
field. The beam is shown in the bottom left corner. Contours are 
228 km s-i (the systemic velocity), ±10, 20, 30, . . . km s'^. 
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3 DATA CUBE MODELLING 

3.1 Velocity field: a first analysis 

From the above-mentioned data cubes at two different res- 
olutions velocity fields were constructed based on the Mod- 
ified Envelope Tracing method (Gentile et al. 2004). A first 
attempt to understand the kinematics of NGC 3741 was 
made by performing a tilted-ring fit (Begeman 1989) to 
the abovementioned velocity fields using the task ROTCUR 
within GIPSY (van der Hulst et al. 1992), and leaving the 
inclination and position angles as free parameters. This way, 
the orientation angles were very poorly constrained. In order 
to obtain a first estimate of the radial run of the position 
angle, we fixed the inclination at the value of 64°, the aver- 
age between the two extremes of the inclination range fitted 
by Begum et al. (2005) ; the results are shown in Fig. |3] The 
kinematical position angles found this way do not agree with 
the morphological position angles, as can be seen in Figs. |4] 
andO In Fig.|4]we compared a model velocity field based on 



Figure 5. Comparison between the observed high-resolution to- 
tal intensity map and the one computed from a model data cube 
based on the output parameters from the tilted-ring fit on the 
velocity field. The beam is shown in the bottom left corner. Con- 
tours are (3, 9, 27) X 10^° atom cm"^. 



the parameters derived from the tilted-ring modelling to the 
observed velocity field: even though the observed and mod- 
elled contours are very close to each other, the overall shape 
of the modelled galaxy does not coincide with the observed 
one. This points towards a misalignment between the mor- 
phological and kinematical position angles and it is a first 
hint at the presence of non-circular motions, that were not 
taken into account by Begum et al. (2005). 



3.2 Improved orientation parameters 

As we have shown in the previous section, the orientation 
parameters derived from the tilted-ring fit of the velocity 
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Figure 6. Comparison between the observed high-resolution 
total-intensity map and that computed from the model data cube, 
built with our final geometrical parameters. The beam is shown in 
the bottom left corner. Contours are (3, 9, 27) X 10^0 atom cm . 
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Figure 7. Comparison between the observed total-intensity map 
of the low-resolution cube and the one computed from the model 
data cube built with our final geometrical parameters. The beam 
is shown in the bottom left corner. Contours are (1.5, 4.5, 13.5, 
...) X 10^^ atom cm~^. 



field cannot account for the observed data cube. We then 
built synthetic data cubes (using a modified version of the 
task GALMOD in GIPSY, Barbieri et al. 2005), which were 
compared to the actual observations. These model observa- 
tions are based on a number of geometrical and kinematical 
parameters that describe the HI disk, such as the rotation 
curve and the inclination and position angles as a function of 
radius. The other input parameters are: the central position 
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Figure 8. Total intensity maps of the low-resolution cube. Ob- 
servations: black contours; Model: grey contours. Top left: using 
our best inclination minus 5°. Top right: our best inclination plus 
5°. Bottom left: our best position angles minus 5°. Bottom right: 
our best position angles plus 5 °. 



and the systemic velocity, which can be derived with reason- 
able accuracy from the observations; the velocity dispersion 
as a function of radius, which we kept at a fixed value of 6 
km s~^, typical of dwarf irregulars (Weldrake et al. 2003); 
the surface density as a function of radius, which is derived 
from the total HI map; the scale height of the HI disk, which 
we assumed to be 300 pc, a typical value for dwarf galaxies; 
the radial motions as a function of radius. With the choice 
of parameters from Section 13.11 the position angle was in- 
correct, and the inclination was underestimated, as can be 
seen in Figs. |4] and (5] These two parameters were then it- 
eratively changed in order to give a better representation 
of the total HI map; the results are show in Figs. |S] and [7] 
for the high- and low-resolution data cubes respectively, in 
which the models clearly match the data much better. Com- 
pared to the position angles shown in Fig. |3] the new ones 
are different by values ranging from 70° in the inner parts 
to about 15° in the outer ones. The inclination was changed 
to a fixed value of 70°, in order to obtain good agreement 
with the total HI maps, both at high and low resolution. 
The radial dependence of the position angle is shown in the 
middle panel of Fig. [121 In Fig. |8] we show low-resolution 
total intensity maps with different choices of the orientation 
parameters. Based upon these models, we estimate the un- 
certainty of the inclination and position angles to be about 
4° . In the innermost few data points we conservatively dou- 
bled this value, given the smaller number of beams that 
sample the total HI map and the higher complexity of the 
HI kinematics and distribution. Note that the assumption 
of a constant velocity dispersion of 6 km s~^ is consistent 
with the observations: even though it is apparently increas- 
ing towards the centre, a velocity dispersion map of a model 
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Figure 9. Median values of tlie velocity dispersion from the ve- 
locity dispersion maps integrated over ellipses. Full circles: from 
the observed high-resolution data cube. Empty circles: from a 
high-resolution model data cube with a constant velocity disper- 
sion of 6 km . 
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Figure 10. Comparison between the observed central channel 
maps of the high-resolution data cube and the modelled ones, us- 
ing the same inclination and position angles as in Fig. |6] and dis- 
regarding non-circular motions. Contours are -2.25, 2.25 (2.5 cr), 
6.75, and 20.25 mjy/beam. The beam is shown in the bottom left 
corner, and the cross indicates the centre of the galaxy. 



data cube with a constant velocity dispersion of 6 km s^^ 
exhibits a very similar behaviour (Fig. |9}, because of the 
broadening of the velocity profiles due to the beam where 
the velocity gradients are the largest. 



3.3 Harmonic decomposition of the velocity field 

In order to account for the kinematics, and therefore to ex- 
plain the difference between morphological and kinematical 
position angles, one must resort to non-circular motions. A 
data cube built with the inclination and position angles that 
match the total HI map is inconsistent with the observa- 
tions if non-circular motions are kept to zero (Fig. llOf) . An 
estimate of their amplitude can be found by performing a 
harmonic decomposition of the velocity field (Schoenmakers 
et al. 1997, Wong, Blitz & Bosma 2004), keeping the po- 



sition and inclination angles fixed at the values derived in 
the previous Section, since leaving them as free parameters 
leads to unconstrained fits. With this method the observed 
hne-of-sight velocity Vios is modelled by 



VloB ~ Co + 



'cjCoa{jTp) 



(1) 



where tp is the azimuthal angle. In this analysis we consid- 
ered terms up to j = 3. The results are illustrated in Figs. 1111 
for the high-resolution cube (the low-resolution cube gives 
similar results), where the si component (indicative of ra- 
dial motions) is non-zero and varies approximately between 
5 and 10 km s^^, with a region between 150" and 250" 
where they are close to 5 km s~^. The very small values of 
C2 and S2 indicate that the kinematics of NGC 3741 is very 
symmetric. Wong et al. (2004) show that an error in the fit- 
ted inclination angle causes the cs term to be proportional to 
ci sin i, which seems to be the case here. However, as seen in 
Fig. [51 the inclination is quite well constrained from the to- 
tal HI map. We interpret this spurious ca oc Ci sini term as 
the difference between an analysis that assumes an infinitely 
thin disk (the harmonic decomposition of the velocity field) 
and one that takes this effect into account (the data cube 
modelling). Indeed, within the framework of the harmonic 
analysis this term would indicate an incorrect value of the 
inclination. However, the harmonic decomposition assumes 
an infinitely thin disk: if the disk thickness is actually non- 
negligible, this assumption will result in an underestimated 
inclination angle. Therefore, if one fixes the inclination at 
its correct value (accounting for the thickness of the disk), 
the output harmonic decomposition procedure will display 
a spurious "ca cx ci sini" effect. 

Note that the plotted errors in Fig. [11] are formal errors 
coming from the fitting procedure and they are likely to be 
smaller than the real uncertainties in the parameters (see 
e.g. Bureau & Carignan 1999). 



3.4 Model data cubes: high resolution 

Starting from the results of the harmonic decomposition of 
the velocity field we built different model data cubes, iter- 
atively improving the input parameters of the models until 
acceptable results were found. Some parameters were ad- 
justed at an earlier stage (the position and inclination an- 
gles), while others (the rotation and radial velocities) needed 
some additional small changes. An automatic and objective 
method to fit a model data cube to an observed one, and 
to estimate the uncertainties in the fit is still a matter of 
active research (Jozsa et al. 2006, in prep.). However, as a 
first estimate of the kinematical properties we used the most 
complete objective method available, the harmonic decom- 
position of the velocity field, whose results were improved 
by matching visually some key features in the model data 
cubes to the observed ones. It is possible to associate a to 
each model but its value cannot be considered for statistical 
purposes, it can only be a useful tool to compare different 
models. The reasons are that 1) there is no objective way to 
define how big the data cube should be or which portion of 
the cube should be fitted and 2) the dumpiness and inhomo- 
geneity of actual distribution of gas cannot be reproduced 
in the smooth axisymmetric models shown here. 
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Figure 11. Results of the harmonic decomposition of the high-resolution velocity field using the geometrical parameters based on data 
cube modelling. 



The rotation curve was derived from the tilted-ring fit 
of the velocity field by only considering those points that are 
less than 45° away from the major axis, while for the radial 
velocity as a function of radius the whole velocity field was 
used. It turns out that the value of the rotation and radial 
velocities had to be slightly modified in the inner parts (com- 
pared to the output of the tilted-ring fit to the velocity field) 
in order to reach a better agreement with the observations. 
The rotation velocity had to be increased by up to 4 km s^^ 
in the innermost parts, and to be decreased by 2 km s~^ in 
the outer parts of the low-resolution cube. The radial veloc- 
ity was increased by up to 4 km s~^ in the innermost and 



outermost parts. The comparison between our best model 
high-resolution data cube and the observed one is shown in 
Fig. 1131 With our choice of input parameters we are able to 
reproduce all the major features present in the observations. 
We reiterate that without the introduction of non-circular 
motions the good agreement would not have been reached, 
due to the impossibility of reproducing simultaneously the 
total intensity map (Figs. [1] and [2]) and the central channels 
(Fig. [To}. The excellent agreement between the two data 
cubes is also visible in Fig. 1141 in which position- velocity di- 
agrams along the major axis are compared. The only feature 
that is not well reproduced is the northern part of the emis- 
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Figure 12. The geometrical and physical parameters that lead to 
the final data cubes. The radial velocities are discussed in Sections 
13.31 13.41 13.51 and |4] The position angle is derived in Section 13.21 
while the rotation curve is derived in Sections 13.31 13.41 and 13.51 



sion in channels 237.5 km to 262.2 km s^^ (receding 
side, see Fig. I13|l . However, since the corresponding emis- 
sion is not present in the approaching channels, it means 
that it is an asymmetric feature of the galaxy, thus it is not 
reproducible with our axisymmetric models. 

Bearing in mind the caveats discussed above on the use 
of the statistics, the fact that our preferred model is a bet- 
ter reproduction of the observations (based on the masked 
cube, see Section [2| can also be seen in the smaller x value 
ix^ = 2.45 X 10* for about 2.1x10* de grees of freedom) com- 
pared to the model data cube built from the first attempt to 
derive the rotation curve (Section 13.11 x^ ~ 3.62 x 10*) or 
compared to the model with the correct orientation angles 
but no non-circular motions (Fig. [TOl = 3.97 x 10*). 

A similar procedure was performed to fit the low- 



resolution cube (see Section [33} ; the results from the low- 
resolution cube were used for radii larger than 180", because 
of its better signal-to-noise ratio in the outer parts. 



3.5 Model data cubes: low resolution 

One of the main goals of this paper is to investigate the 
distribution of dark matter at large radii; we smoothed the 
high-resolution cube to obtain a beam of 60" x 60" in order 
to highlight the faint diffuse emission. 

Similarly to the high resolution cube, first a tilted-ring 
fit on the velocity field was performed. Again, the kine- 
matical position angle is different from the morphological 
one, and it needs to be modified in order to reproduce the 
data cube. The same procedure as above was applied to 
look for the best geometrical and physical parameters to 
describe the observed data cube. The parameters for the 
inner parts were fixed at the values derived from the high- 
resolution cube; based on the signal-to-noise ratio of the 
high-resolution cube, it was decided to only use the high- 
resolution values up to 180". A harmonic decomposition of 
the velocity field was performed, whose results are not shown 
because they are similar to Fig. 1111 The Si term in the outer 
parts is close to 10 km s^^, while the non-zero C3 term is 
again interpreted as the effect of disk thickness; the inclina- 
tion is well constrained from the total HI map (Fig. [T}. 

The comparison between the observed and modelled 
low-resolution total HI maps is displayed in Fig. [T] Again, 
the agreement between the maps is excellent. As can be 
seen in Figs. [TS] and \Wi our choice of parameters enables 
us to reproduce the observations in detail. Without the in- 
troduction of non-circular motions such a good agreement 
would not have been possible. The final choice of inclina- 
tion, position angle, rotation velocity and radial velocity is 
shown in Fig. [121 The parameters from the low-resolution 
cube were sampled every 30" (half the beam FWHM), so 
that the low- and high-resolution cubes contribute to the 
final rotation curve with a comparable number of points. 
Thanks to the exceptional extension of the HI disk, the ro- 
tation curve extends out to 6.6 kpc, which correspond to 42 
B-band exponential scale lengths (see Fig. [2} . 

Given the small size of the optical disk (R25 = 
30'.'7, taken from LEDA, the Lyon-Meudon Extragalactic 
Database) compared to the beam, and given that the ori- 
entation of the optical disk is not very well defined (e.g. 
the position angle estimates of the optical disk range from 
-15.5°, LEDA, to -1-23°, Vaduvescu et al. 2005), it is hard 
to define the radius at which the warp commences. Accord- 
ing to Briggs (1990), warps usually start around R25. For 
the final rotation curve we used the values derived from the 
high- resolution cube out to 180", beyond which the low- 
resolution cube was used. 



3.6 Asymmetric drift 

We applied the following correction for asymmetric drift, in 
order to derive the circular velocity Vc from the observed 
rotation velocity Kot: 



9 In p d \n 



9 In r 9 In r 



(2) 
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Figure 13. Comparison between the observed high-resolution data cube (black contours) and the modelled one (grey contours). Every 
second channel map is shown. Contours are -2.25, 2.25 (2.5 a) and 6.75 mjy/beam. The beam is shown in the bottom left corner, and 
the cross indicates the centre of the galaxy. 
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Figure 14. Comparison between position-velocity diagrams par- 
allel to the major axis of the observed high-resolution cube (black 
contours) and the modelled ones (grey contours). The slices 
are spaced by 70". Contours are -1.8, 1.8 (2 a), 4.5, and 11.2 
mjy/beam. 



where a is the one-dimensional velocity dispersion of 
the gas and p is the gas density, which can be derived from 
the observed surface density distribution under the assump- 
tion of a constant scale height. As shown in Section 13.21 
the observed velocity dispersion in NGC 3741 is consistent 
with a constant value of 6 km s~^ , which is what we assumed 
here. The asymmetric drift correction turns out to be smaller 
than the error bars except for the outermost two data points, 
where the corrections are respectively 2.7 and 4.8 km s~^. 
Note that this correction is very uncertain, mainly because 
we are ignoring the other gaseous components of the ISM, 
in particular the molecular gas. 

In order to estimate the uncertainties in the rotation 
curve, tilted-ring fits to the velocity fields were made such 
as those illustrated in Figs. 1111 but considering one side of 
the velocity field at a time. The errors of the final rota- 
tion curve were taken to be the geometric sum of half the 
difference between the approaching and receding sides and 
half the asymmetric drift correction. A minimum error of 
2/sini km s~^ (half the velocity resolution) was assumed. 



4 NON-CIRCULAR MOTIONS 

The distribution and kinematics of the HI in NGC 3741 
could only be modelled by resorting to non-circular motions 
(Sect.[3|. This is therefore a caveat in the analysis of the dark 
matter distribution, in that the gas is not in perfect circular 
motion and the radial component of its velocity, especially 



in the inner parts, is non-negligible. In the inner ~30" the 
non-circular motions that we found are of the same order 
as the rotation velocity. The innermost points of the rota- 
tion curve are also those that give the strongest constraints 
on the steUar mass-to-light ratio; therefore, the best-fit val- 
ues inferred in the next section are have larger uncertainties 
than the formal ones, since the rotation velocity Kot is not 
likely to be a fair tracer of the circular velocity Vc ■ The non- 
circular motions as derived from the harmonic decomposi- 
tion are around 10 km s~^ between the centre and 80", then 
they drop to about 5 km s~^, between 150" and 200" they 
increase again around 10 km s~^. As discussed in Sect. |3l 
the higher values were increased by a few km s~^ in order 
to better reproduce the observed data cube, so the range of 
radial velocities becomes 5. . .13 km s~^. Some dwarf galax- 
ies (such as NGC 4605, Simon et al. 2005) exhibit a similar 
behaviour in terms of non-circular motions, while in others 
(e.g. DD047, Gentile et al. 2005) their amplitude is much 
smaller. However, the interpretation of these motions is not 
obvious. 

A rough idea can be obtained from the parameters si 
and S3 of the harmonic decomposition (Wong et al. 2004). 
In the inner regions there is a hint of an anti-correlation 
between the two parameters, which can be an indication of 
elliptical streaming. The most likely cause for this kind of 
motion is an inner bar; Figs. [1] and [6] also suggest the pres- 
ence of a bar at the centre. In the outer parts the S3 term 
is very close to zero, which is consistent with a radial fiow. 
In order to disentangle between outfiow and infiow we need 
to know the true orientation of the galaxy, i.e. whether it is 
rotating clockwise or counterclockwise. In the map of total 
HI intensity (Fig. [T} there is a hint of spiral structure. As- 
suming that it is really present and that the spiral arms are 
trailing, we can infer that the galaxy is rotating clockwise, 
so the radial fiow could be interpreted as an infiow. To try 
to better understand these radial motions, we separated the 
galaxy in two parts, dividing along the major axis (because 
the effect of radial motions is stronger along the minor axis). 
The radial motions on the north-western side are higher (by 
~7 km s~^) than on the opposite side, which makes the ac- 
cretion hypothesis more likely. However, that the magnitude 
of the motions seems rather high compared to the theoret- 
ical expectations (e.g. Blitz 1996 expects an infiow velocity 
of 7 km s~^ in the outer HI disk of the Milky way). We 
also note that a full investigation of this phenomenon goes 
beyond the scope of this paper. 

With this assumption, it might be that the material in 
the outermost regions that is seen having a ~ 10 km s~^ 
infiow is being accreted onto the galaxy. The consistency of 
the HI observations with an inner bar and outer accretion 
can be shown by comparing the Si and S3 terms; Wong et 
al. (2004) showed that they can be use to distinguish be- 
tween a bar (elliptical streaming), where the two terms are 
anticorrelated, and a radial fiow, where Si is different from 
zero and S3 is close to zero. We plot the two terms (nor- 
malised to ci) in Fig. 1171 separately for the inner regions 
(radius smaller than 180", using the high-resolution data) 
and for the outer regions (radius larger than 180", using the 
low-resolution data). In the inner parts the two terms are 
anticorrelated, suggesting the presence of a bar, while in the 
outer parts they are uncorrelated, which is consistent with 
a radial (in)fiow. 
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Figure 15. Comparison between the observed low resolution observed data cube (black contours) and the modelled one (grey contours). 
Every second channel map is sliown. Contours are -3.25, 3.25 (2.5 a), 13, and 52 mjy/beam. The beam is shown in the bottom left 
corner, and the cross indicates the centre of the galaxy. 
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Figure 16. Comparison between tlie major-axis position-velocity 
diagrams of the observed low-resolution cube (black contours) and 
the modelled one (grey contours). The slices are spaced by 70". 
Contours are -3.25, 3.25 (2.5 cr), 13, and 52 mjy/beam. 



5 MASS MODELS 

The data presented in the present paper allow us to make 
use of the most extended rotation curve of a galaxy ever 
measured, in terms of the optical disk scale length. The ob- 
served rotation curve Vohs, corrected for asymmetric drift, 
was decomposed into the stellar disk, gaseous disk and dark 
halo contributions Vdisk, Vgas and Vhaio via 



(3) 



The shape of the stellar disk contribution is derived 
from the NIR photometry presented in Vaduvescu et al. 
(2005). From the K-band photometric profile the contribu- 
tion of the stellar disk was derived using the task ROTMOD 
in GIPSY with the assumption of a vertical sech^ distribu- 
tion having a scale height zo — h/6 (van der Kruit & Searle 
1981). Its amplitude is scaled using a constant mass-to-light 
(M/L) ratio. From the photometric profile we can exclude 
any significant contribution from a bulge, because NGC 3741 
even shows a brightness deficit at small radii, extrapolating 
the exponential fit to the inner region. The gaseous disk con- 
tribution was derived from the HI observations presented 
here, multiplying the neutral hydrogen surface density by 
1.33 to account for primordial He. 



5.1 Stellar disk 

The stellar A'l/L ratio being one of the largest uncertainties 
in the rotation curve fits, we describe here the attempts we 
made to put some strong a-priori constraints. Unfortunately, 
with the present data this was not possible to achieve. At- 
tempts were made using the measured (B-K) colour and the 
predictions of stellar population synthesis models of Bell & 
de Jong (2001). We find a predicted K-band stellar M/L 
ratio of 0.3. We also tried to constrain the stellar M/L ra- 
tio by using the predictions of Bell et al. (2003). They only 
give values based on the observed B-V and B-R colours. By 
taking the latter value from LEDA (Lyon-Meudon Extra- 
galactic Database) we obtain a K-band stellar M/L ratio of 
0.7. 

Another attempt at constraining the stellar M/L ra- 
tio was made through an estimate of the lower limit of the 
stellar mass by comparing UBV photometry (NED) plus JK 
magnitudes (Vaduvescu et al. 2005) with several libraries of 
SEDs computed with the spectral synthesis code GRASIL 
(Silva et al. 1998). The observed SED can be well repro- 
duced by several difi'erent star formation histories. There 
are no difficulties in producing an amount of about 10^ M© 
or more in stars as predicted in the cored and in the NFW 
mass models. Instead, in order to obtain significantly lower 
masses, i.e. 5 x 10^ Mq, and a good fit to the observed 
SED, there are some requirements to be fulfilled: a peculiar 
SFR, a delayed formation for 10 Gyr and most of the stars 
formed in the last 0.5 Gyr. In any case, also accepting this 
fact, by adopting a standard Salpeter IMF extended from 
0.1 to 100 Mq, we never found that the observed SED im- 
plies stellar masses lower than ^ 5 x 10^ Mq. These results, 
combined with the estimates from the stellar population syn- 
thesis models, induced us to consider in the rotation curve 
fits a minimum K-band AI/L ratio of 0.3, corresponding to 
a stellar mass of ~ 7 x 10^ Mq at a distance of 3 Mpc. 

We remind, however, that the stellar M/L ratio com- 
ing out of the fits is mostly constrained by the innermost 
two/three data points of the rotation curve, which are also 
the most affected by non-circular motions, in proportion to 
the rotation velocity. In the specific example shown in Fig. 
[18] the stellar mass is ~ 6 x 10^ Mq. The spectrophotomet- 
ric determination of the stellar disk mass does not help in 
discriminating between cored and NFW halos, and it is not 
decisive against MOND. In any case, the small stellar mass 
found with these estimates (of the order of lO^M©) is an 
advantage for our aims: uncertainties in the stellar M/L do 
not have a great infiuence on the derived dark matter halo 
profile. 



5.2 Dark halo 

Concerning the dark matter halo, several alternatives were 
investigated. In numerous previous studies dark halos with 
a central constant-density core provided the best fits to the 
rotation curve. An example of such a cored halo is the Burk- 
ert halo (Burkert 1995, Salucci & Burkert 2001); its density 
distribution is given by 



PBur{r) 



3 

PO^corc 



^corc ) 



(4) 
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were po (the central density) and rcorc (the core radius) are 
the two free parameters. 

Then, another fit to the observed rotation curve was 
performed using the Navarro, Frenk and White halo (NFW) , 
the result of an analytical fit to the density distribution of 
dark matter halos in CDM cosmological simulations. 



Pnpw(»") 



(r/r,)(l + r/r,)2' 



(5) 



Here, pa and r-g are the characteristic density and scale of the 
halo. More recent simulations (Navarro et al. 2004) have im- 
proved the NFW fit to the ACDM halo density distribution, 
but on scales of radii we are considering here the NFW halo 
is still a good fit to the outcome of the ACDM simulations. 
The two parameters ps and rs are in principle independent, 
but they were shown to be correlated in simulations (e.g. 
Wechsler et al. 2002) through the virial mass A/vir, the con- 
centration parameter Cvir = rvir/ra and the critical density 
of the Universe pcrit- The recent third year WMAP results 
(Spergel et al. 2006) suggest however that the concentration 
parameters predicted by Bullock et al. (2001) and Wechsler 
et al. (2002) should be lower, in particular due to the lower 
as. As the "original" Cvir — Mvir relation we considered eq. 
18 in Bullock et al. (2001), whose normalisation was lowered 
by 25%, compatible with 20% used in the model of Button 
et al. (2006), 28% that comes out of Bullock et al.'s CVIR 
code and 25% found by Gnedin et al. (2006). The relations 
we used were therefore: 



13.6 



lOiiM 



8.8 



M ■ 
IQiiMc 



r Pcrit 



kpc(6) 



(7) 



3 ln(l + Cvir) — Cvir/(1 + Cvir) 

where A is the virial overdensity at z = that can be 
computed following Bryan & Norman (1998). The scatter 
in the relation that links Cvir to Mvir has been reported to 
be Alogc ~ 0.14 (Bullock et al. 2001), although this value 
slightly varies from one study to another. We also note that 
these relations assume a running of the tilt of the power 
spectrum {dus / dink) of 0, while different values lead to sig- 
nificant changes in the normalisation of these relations. 

5.3 MOND 

A hypothesis that has been put forward (Modified Newto- 
nian Dynamics, MOND; Milgrom 1983) to explain the ab- 
sence of a Keplerian decline in rotation curves is that New- 
tonian gravity does not hold below a certain acceleration ao, 
without the need of invoking dark matter. The true gravi- 
tational acceleration g and the Newtonian one qn are linked 
through: 



9 ' 



9n 



pig/ao) 



(8) 



where I-l{x) is an interpolation function that tends to 1 
for g ^ ao and tends to g/ao for g <^ ao. Traditionally fi{x) 
was given the following form: 



(9) 



Recently, Famaey & Binney (2005) proposed the follow- 
ing alternative form of ij.{x), which, contrary to Eqn.|9l Zhao 



& Famaey (2006) have shown not to be in contrast with the 
relativistic MOND theory of Bekenstein (2004): 



li{x) = 



1 + 2: 



(10) 



In both cases, the stellar M/L ratio and the distance 
were left as free parameters. 



5.4 Gas scaling 

A noticeable global property of rotation curves of disk galax- 
ies is that the ratio of the gaseous disk and dark matter sur- 
face density is approximately constant with radius (Hoekstra 
et al. 2001). This implies that one can attempt to fit a rota- 
tion curve without a dark matter halo by simply scaling up 
the contribution of the gaseous disk. This fact is important 
in view of the hypothesis that a substantial fraction of dark 
matter within the radii reached by rotation curves might 
reside in a disk, in the form of cold clumps of molecular 
hydrogen (Pfenniger, Combes & Martinet 1994). Previous 
studies have shown that this hypothesis does not reproduce 
the details of observed high-quality rotation curves of spiral 
galaxies (Gentile et al. 2004), even though it is a remarkable 
feature of disk galaxies. 



6 RESULTS 

The results of the mass decompositions are shown in Fig. [19] 
and in Table 1. The Burkert halo gives a very good-quality 
fit, with a core radius about five times larger than the op- 
tical radius. Its best-fit AI/L ratio is quite high, rather 
poorly constrained and consistent with the spectrophoto- 
metric analysis based on the GRASIL code. Also, if the 
M/L in the B-band is considered, we find very similar values 
compared to that found by Carignan & Purton (1998) for 
DDO 154, whose properties (e.g. very large extent of the HI 
disk, maximum rotation velocity) are similar to NGC3741. 
To further test whether reasonable stellar M/L may weaken 
the support for a cored halo we made a Burkert fit with a 
fixed M/Lk = 1 (within the range derived by Bell et al. 
2003 for blue galaxies): we find that the fit is still very good 
(xL = 0.62). 

The NFW fit using Eqn.|S]is poor, with the usual effect 
of overestimating the dark mass density in the inner parts 
and underestimating it in the outer parts. Moreover, the 
best-fit stellar M/L ratio is 0: more realistic values of M/L 
would only make the fit worse. 

If the relation between Cvir and Afvir is relaxed, a good 
fit can be made. However, this is achieved at the expense 
of rather unrealistic (in the context of ACDM) best-fit pa- 
rameters: the best-fit Cvir is 6, which is about 2.5 a off the 
predicted relation (corrected for the third year WMAP cos- 
mological parameters). This means that the rotation curve 
of NGC 3741 can be fitted with an NFW halo but only if it 
is a 2.5-cr exception to the Cvir- A/vir relation. Also, note that 
in this case the NFW best-fit virial mass is too high (~ 10^^ 
M0) for a galaxy with an observed maximum velocity of 
45. . .50 km s~^, even though the constraint on M^ir is weak. 
Moreover, a Afstars /Afvir ~ 10"'' is difficult to explain even 
taking the feedback from SN explosions into account. 
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If the distance is kept as an unconstrained free parame- 
ter, the standard MOND fit is almost as good as the Burkert 
fit, and so is the MOND fit with the new ^(a;) (eq. [TO]). In 
both MOND fits the acceleration parameter ao was kept 
fixed at i.2 X fO~* cm s~^ (Begeman, Broeils & Sanders 
199f). The best-fit distance in the MOND fits is significantly 
higher than the value quoted by Karachentsev et al. (2004) , 
3.0 ± 0.3 Mpc, but it is consistent with the distance of 3.5 ± 
0.7 Mpc given in Georgiev et al. (f997). The former is more 
likely to be a reliable distance estimate, as it was made us- 
ing the tip of the Red-Giant Branch instead of the brightest 
blue stars, so the high best-fit distance might potentially 
constitute a problem for the standard MOND jJ.{x) (see Ta- 
ble [TJ. By allowing the distance to vary only within the 3.0 
± 0.3 Mpc range we found that the standard ^i{x) gives a 
relatively poor fit (x?od = 1.41, which is however still much 
better than the constrained NFW fit), while the new fi{x) 
gives a better fit (x?od = 0.91); this fit can be further im- 
proved (x?ed = 0.64) if ao is assumed to be 1.35xl0~* cm 
s~^ (see Famaey et al. 2006). 

Scaling up the contribution of the gaseous disk gives a 
reasonably good fit, with a scaling factor for the gas surface 
density of ~ 13, in agreement with the results of Hoekstra 
et al. (2001). 

A fit of the rotation curve was also made by assuming 
that the dark matter density can be described as a power-law 
(e.g. Simon et al. 2003), defined as p oc r". This is not the 
case in most models of dark matter halos; a power-law fit will 
therefore provide a sort of average slope of the dark matter 
density profile over the range of radii probed by the rotation 
curve. Note that the outcome of such a fit depends on the 
extension of the rotation curve: a more extended rotation 
curve will have a steeper average slope. We fixed the stellar 
mass-to-light ratio at the average value between the two 
extremes found in the above fits: M/L = 0.9. This quantity 
turns out not to have a great influence on a, since the disk 
dominates the kinematics only in the innermost few points. 
We find Q = —1.11 ± 0.11 and Xred ~ 0.35, which however 
does not mean that NGC 3741 has a cuspy dark matter 
halo. Indeed, at their respective core radius, a Burkert halo 
has a dlogp/dlogr — —1.5, and a pseudo-isothermal halo has 
dlogp/dlogr = —1.0. This shows that if the rotation curve 
extends beyond the core radius (such as in NGC 3741), the 
average slope, as measured by a power-law fit, will be much 
steeper than zero. On the other hand, for an NFW halo it 
is difficult to reproduce an average slope close to -1, since 
its slope is already -2 at r = (which is of the same order 
as the last measured radius of our NGC 3741 observations). 
See e.g. Fig. 3 of Navarro et al. (2004) for the variation of 
the density slope with radius in ACDM simulations of dwarf 
galaxies. The only ways to reconcile an NFW profile with the 
observations are to have a very large Vs, which means that 
either the virial mass has to be very high, or to decrease the 
concentration, which is at odds with the ACDM predictions 
in the concordance cosmology. 

NGC 3741 has therefore evidence for a cored halo. This 
could be linked to the presence of a bar. Indeed, barred 
galaxies seem to have somewhat shallower dark matter den- 
sity profiles (Swaters et al. 2003), despite the uncertainties 
in the mass models due to non-circular motions. However, 
there is no general agreement yet as to whether realistic bars 



Table 1. Best-fitting results: parameters and associated la uncer- 
tainties, reduced ■ ''core is in kpc, po is in units of 10~^* g cm"^, 
Mvir is in units of M©, and the distance is in Mpc. Mstars and 
Mgas are the masses of the stellar and gaseous disks, respectively, 
given in Mq. Afi,ar=A^stars-l- Afgas is the baryonic mass. Afvir in 
the case of the Burkert halo was computed by extrapolating the 
density profile until the average density reached A Pcrit i where A 
is 92.5. 
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can significantly flatten the cusp (Weinberg & Katz 2002, 
Sellwood 2006)! 



7 CONCLUSIONS 

We have presented the analysis of HI observations of the 
nearby dwarf irregular galaxy NGC 3741, performed at the 
WSRT. The HI disk of this galaxy is very extended, the most 
extended ever observed in terms of optical scale length. The 
huge extension of the HI disk enables us to trace the rotation 
curve out to unprecedented distances in terms of the size of 
the optical disk: the last radius at which we can trace the 
kinematics is 42 B-band exponential scale lengths, approx- 
imately 15 times R25. The HI disk displays a warp, which 
is very symmetric. In order to derive the distribution and 
kinematics, we built model data cubes; with our choice of 
geometrical and kinematical parameters, the observed data 
cube is reproduced accurately. Some key features in the data 
cube are characteristic of non-circular motions, whose am- 
plitude turns out to be of the order of 5. . .13 km s~^. Their 
interpretation is non-trivial, but they are consistent with an 
inner bar and accreting outer material. 

Subsequently, the standard rotation curve decomposi- 
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Figure 17. Comparison between the si and S3 parameters of the harmonic decompositions of the velocity fields. Left: for the high- 
resolution data and radii smaller than 180". Right; for the low-resolution data and radii larger than 180". 
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Figure 18. Observed (stars) and fitted (solid line) SEDs, using 
the GRASIL code and a stellar mass of 6 X lO^Mg. 



tion was performed. The (cored) Burkert dark halo fits the 

data very well; the fit performed by taking the dark matter 
density distribution predicted by the ACDM cosmological 
simulations fits badly, unless the prediction of a correlation 
between the concentration parameter and the virial mass is 
relaxed. In this case, however, the price to pay is a concen- 
tration parameter 2.5 a below the predicted c— Mvir relation 
and a high virial mass (but poorly constrained) of 10^^ Mq. 

Scaling up the contribution of the gaseous disk also 
gives acceptable fits. MOND, both with the standard inter- 
polation function and with the new one proposed by Famaoy 
& Binney (2005), is also consistent with the observed rota- 
tion curve. We also notice that the distance to NGC 3741 
might be a potential problem for the standard MOND in- 
terpolation function, but not for the new one. 
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